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ABSTRACT 

BL Lacertae (Lac) objects that are detected at very-high energies (VHE) are of fundamen- 
tal importance to study multiple astrophysical processes, including the physics of jets, the 
properties of the extragalactic background light and the strength of the intergalactic mag- 
netic field. Unfortunately, since most blazars have featureless optical spectra that preclude a 
redshift determination, a substantial fraction of these VHE extragalactic sources cannot be 
used for cosmological studies. To assess whether molecular lines are a viable way to estab- 
lish distances, we have undertaken a pilot program at the IRAM 30m telescope to search for 
CO lines in three BL Lac objects with known redshifts. We report a positive detection of 
Mh 2 ~ 3 x 10 s M Q toward 1ES 1959+650, but due to the poor quality of the baseline, this 
value is affected by a large systematic uncertainty. For the remaining two sources, W Comae 
and RGB J0710+591, we derive 3cr upper limits at, respectively, Mh 2 < 8.0 x 10 s M Q and 
Mh 2 < 1.6 x 10 9 M Q , assuming a line width of 150 km s _1 and a standard conversion factor 
a — 4 M (Kkms -1 pc 2 )^ 1 . If these low molecular gas masses are typical for blazars, blind 
redshift searches in molecular lines are currently unfeasible. However, deep observations are 
still a promising way to obtain precise redshifts for sources whose approximate distances are 
known via indirect methods. Our observations further reveal a deficiency of molecular gas in 
BL Lac objects compared to quasars, suggesting that the host galaxies of these two types of 
active galactic nuclei (AGN) are not drawn from the same parent population. Future obser- 
vations are needed to assess whether this discrepancy is statistically significant, but our pilot 
program shows how studies of the interstellar medium in AGN can provide key information 
to explore the connection between the active nuclei and the host galaxies. 

Key words: BL Lacertae objects: general, BL Lacertae objects: individual: 1ES 1959+650, 
galaxies: ISM, galaxies: distances and redshifts, galaxies: active, ISM: molecules. 



1 INTRODUCTION 

Blazars are among the most violently variable radio sources within 
the population of active galactic n uclei (AGN). In the classical uni- 
fied model (e.g. lAntonuccil 1 19930 . they are considered the face- 
on view of radio-loud AGN, with a beamed jet pointing along 
the observer line of sight and with blobs of highly relativistic 
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particles moving down the jet. The non-thermal radiation pro- 
duces a double-peaked spectral energy distribution (SED), with 
the lower energy peak resulting from the synchrotron radiation 
of relativistic leptons in the presence of a tangled magnetic field, 
while the higher energy peak is commonly attributed to inverse- 
Compton up-scattering by the relativistic particles within the jet 
of either the synchrotron photons themselves (namely synchrotron 
self-Compton emission; SSC) or a photon field external to the jet 
(namely external Compton e mission; EC). A revi e w of these lep- 



tonic models can be found mlpermer et al.l 
dl992h : iMarscher et all fcoOOl) ; ISikora et al 
ences therein. 
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Within this general class, blazars can be further sub-classified 
as flat-spectrum radio-loud quasars (FSRQs) if broad emission 
lines are visible or, otherwise, as BL Lacertae (Lac) objects. BL 
Lac objects are then categorized based on the frequency location of 
the lower energy peak (f sy nch), with low-synchrotron-peaked BL 
Lac (LSP) objects having a ^ sync h below 10 14 Hz, intermediate- 
synchrotron-peaked BL Lac (ISP) objects peaking between 10 14 
and 10 15 Hz and high-sychrotro n-peaked BL Lac (HSP) objects 
showing a peak above 10 15 Hz jAbdo et alJlioiObh . Across these 
differ ent sub-classes, bla zars exhibit a continuous variation in their 
SEP jFossati et al.l 1998h and this has been i nterpreted as an evolu- 
tionary sequence (Bot tcher & Der mer 2002) associated to the vari- 
ation of the diffuse radiation field in the surroundings of the rel- 
ativistic jet dGhisellini et alj [l998l) . The evolution along the se- 
quence is connected to the decrease in the accretion r ate and to 
the onset of the r adiatively inefficient accretion flow JHoI 12008k 
iTrump et al.ll201 lh . At the beginning of the sequence, blazars have 
appreciable external radiation field that facilitates effective cooling, 
resulting in low frequency peaked synchrotron and inverse Comp- 
ton emission. A progressively lower contribution of the external 
field along the sequence inhibits efficient cooling, resulting in an 
higher frequency peak of the synchrotron and inverse Compton 
emission. 

New highly sensitive 7-ray instruments such as the VERI- 
TAS and HESS Cherenkov arrays, the Fermi Large Area Telescope 
(LAT), and the MAGIC Telescopes are offering a new perspective 
of the complete electromagnetic spectrum of blazars, providing the 
oppo rtunity to investi gate the origin of the highest energy photons 
(e.g. lAliu etaT]|20T ]]), and to use these powerful AGN as back- 
ground sources to prob e the properties of the extragalactic back- 
ground light (EBL; e.g. lAharonian et al1l2006l) . High energy 7-ray 
photons interact with low energy EBL photons via pair production, 
suppressing the intensity of the 7-ray spectra emitted by extragalac- 
tic objects. Therefore, the absorbed spectrum contains information 
on the EBL that, in turn, encodes the integrated history of struc- 
ture formation and the evolution of stars and galaxies in the Uni- 
verse. Further, these pair-produced electrons lose energy by inverse 
Compton scattering against the cosmic microwave background ra- 
diation. Since the resulting flux is dependent on the intergalactic 
magnetic field (I GMF), 7-rays emitting blaz ars are useful probes 
of the IGMF (e.g. N eronov & Semiko zl l2009h . However, since the 
non-thermal emission from BL Lac objects, by definition, outshines 
the weaker optical spectral features that are suitable for establish- 
ing redshifts, the distances to blazars that are detected at very high 
energies (VHE; E > 100 GeV) are sometimes unkno wn, limiting 
the u se of these objects for cosmological studies (e.g. lAbdo et all 
[20Tl|). 

A possibly powerful yet unexplored way to establish the dis- 
tance to BL Lac objects is the search for bright carbon monox- 
ide emission lines, such as 12 CO(l - 0) and 12 CO(2 - 1). In- 
dee d, previous CO ob s ervations of quasars and Seyfert I galax- 
ies jEvans et alj|200ll ; IScoville et aU 120031 : iBertram et ail i2007l) 
have yielded a high detection rate (~ 70%) of molecular gas 
from the AGN host galaxies. The rat e of detection is lower in 
radio-bright galaxies (~ 60 — 40%; Ocana Flaquer et al. l2010l : 
Smolcic & Riechers 2011), but one in two or three sources is de- 
tected in moderately deep observations. Inspired by these previ- 
ous results, we investigate whether molecular emission lines are in 
fact a promising way to measure distances in BL Lac objects. Be- 
sides this practical reason, establishing the molecular content in BL 
Lac host galaxies compared to other classes of AGN potentially of- 
fers additional insight into the origin of the blazar evolutionary se- 



quence, the unified model for AGN and the processes of feedback 
in massive galaxies. 

This paper presents results from pilot observations of three 
BL Lac objects at the IRAM 30m telescope. The target selection 
for this study is presented in Section [2] while observations and 
the data reduction are discussed in Section [3] Section [4] offers an 
overview of the molecular gas properties of the three blazars tar- 
geted by our observations, while the discussion follows in Section 
[5] The summary and future prospects conclude this paper in Sec- 
tion^ Througho ut this work we ad opt cosmological parameters 
from WMAP7 dKomatsu et aT]|201 ll) . 



2 TARGET SELECTION 

For this pilot study, we searched the literature for VHE detected 
BL Lac objects with a reliable redshift determination and a lu- 
minosity distance Dl < 600 Mpc. Our final selection includes 
RGB J0710+591, W Comae, and 1ES 1959+650, three blazars that 
fall within the class of radio and X-ray bright source s. Details on 
the VHE properties of these sou rces can be found inlAcciari et all 
( l2009h . lAcciari et al.l JioToh . and lAharonian et al.ld2003h . 

RGB J0710+591 is an HSP object at z = 0.125, whose red- 
shift is secured by the detecti on of Ca H & K and G band abso rption 
lines in the optical spectrum dLaurent-Muehleisen et al . 1998). The 
blazar host galaxy is fully res olved and the light profile is charac- 
teristic of an elliptical galaxy dScarpa et al1l2000l) . 1ES 1959+650 



is again an HSP object at z = 0.047 who se redshift is provid ed by 
prominent Ca H & K absorption lines dSchachter etal.ll 1 993). The 
host galaxy is fully resolved and the radial profile is very well de- 
scribed by a point s ource combined with a de Vaucouleurs model 
dScarpa et alJkOOCh . Finally, W Comae is an ISP object, whose 
redshift at z = 0.103 can be established from the emission lines 
[O n], [O in], H a, and [N 11] detected in the Sloan Digital Sky 
Survey spectrum ( Abazaiian et al. 2009). As visible from Figure 
[T] a companion galaxy is located in the south-west direction from 
this source, with the two objects possibly in interact ion. The light 
profil e is well described by a Sersic index n ~ 12 dNilsson et al.l 
120031) . 



3 OBSERVATIONS AND DATA REDUCTION 

We observed the three blazars with the IRAM 30m telescope on 
Pico Veleta (Spain) between UT 15 and 19 July 2011. Additional 
observations for RGB J0710+591 were taken between UT 2 and 
5 September 2011. Weather conditions were generally good, with 
system temperatures ranging between 100 and 150 K at ~ 110 
GHz. 

To acquire the 12 CO(l — 0) line (rest-frame frequency va = 
115.2712018 GHz) and the 12 CO(2 - 1) line (v = 230.5380000 
GHz) simultaneously, we used four single pixel heterodyne EMIR 
receivers, two centered on the E0 band (3 mm) and two on the 
E2 band (1 mm). The local oscillator was tuned such that the red- 
shifted 12 CO(l — 0) lines would fall at the center of the lower inner 
sideband for RGB J0710+591 and W Comae and, due to its lower 
redshift, in the upper inner sideband for 1ES 1959+650. The tele- 
scope half-power beam widths are ~ 23" and ~ 12", respectively 
at the 12 CO(l - 0) and 12 CO(2 - 1) redshifted frequencies, and 
encompass the host galaxies (Figure [T). The data were recorded 
using the WILMA autocorrelator providing a spectral resolution 
of 2 MHz. This corresponds to a resolution of ~ 6 km s _1 at 
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Figure 1. A"-band images from the Two Micron All Sky Survey JSkrutskie et"aT]|2006l) for RGB J0710+591 (left), 1ES 1959+650 (center), and W Comae 
(right). The two dashed and dotted circles around each blazar represent the beam size at the observed frequencies for the 12 CO(l — 0) and 12 CO(2 — 1) 
transitions, respectively. In all cases, the beam encompasses the blazar host galaxies. 



Table 1. Summary of the molecular emission properties of RGB J0710+591, 1ES 1959+650, and W Comae. 









12 CO(l - 0) 


12 CO(2 - 1) 






Name 






7 3 T 4 
^co - 1 rms 


r 5 


rp 4 

-t rms 


T 5 

1 CO 


T 6 
^CO 


M H2 7 




(Mpc) 




(mK/mJy) 


(mjy km s _1 ) 


(mK/mJy) 


(mJy km s — 1 ) 


(10 7 Kkms" 1 pc 2 ) 


(1O 8 M ) 


RGB J0710+591 


583 


0.125 


0.20/1.2 


< 184 


0.44/3.3 


< 491 


< 13.66 


< 5.47 


1ES 1959+650 


208 


0.047 


0.0474 0.18/1.1 


754 ± 170* 


0.44/3.4 


1647 ±511 


7.60 ±1.71* 


3.04 ±0.69* 


W Comae 


473 


0.103 


0.15/0.9 


< 134 


0.30/2.2 


< 334 


< 6.67 


< 2.67 



luminosity distance. 2 Optical redshift. 3 Molecular line redshift. 4 Root-mean-square antenna temperature and flux density at the final resolution of ~ 45 
km s _1 per channel. 5 Line intensity, where for non-detections we report the la limit assuming a line width of 150 km s — 1 . 6 Line luminosity of the 
12 CO(l — 0) transition. 7 Molecular gas mass assuming a = 4Mq (K km s — 1 pc 2 ) - 1 . * The line flux and luminosity for the 12 CO(l — 0) transition are 
affected by the poor quality of the baseline (see Section[3]for details). 



3 mm and ~ 3 km s _1 at 1 mm. The observations were con- 
ducted in wobbler-switching mode, with a switching frequency of 
0.5 Hz and a symmetrical azimuthal wobbler throw of 50" to maxi- 
mize the baseline stability. For each blazar, we acquired series of 12 
ON/OFF subscans of 30 seconds each. Calibrations were repeated 
every 6 minutes, while focus and pointing were checked throughout 
the observing period with bright sources. At the end of the observ- 
ing run, we had acquired 4.8 hours on source for 1ES 1959+650, 
7.6 hours for W Comae and 5.4 hours for RGB J0710+591. 

The data reduction was completed with the IRAM software 
CLASS/GILDAS. The flux calibrated scans were subsequently 
combined and analyzed using a set of in-house IDL procedures. 
To compensate for the variability of the continuum level that is 
common for blazars, we subtracted a constant baseline from each 
scan before coadding. After applying the heliocentric correction to 
homogenize observations at different epochs, we combined all the 
scans, weighting by the exposure time and the system temperature. 
Both polarizations were included in the final stack. 

Due to the bright continuum level at 3 mm for 1ES 1959+650 
and W Comae (~ 0.16 Jy and ~ 0.23 Jy, respectively), the qual- 
ity of the baseline is degraded by the presence of standing waves 
between the secondary and the feed to the receiver. Specifically, 
the spectrum of 1ES 1959+650 at 3 mm exhibits ripples with 
a peak-to-peak amplitude of ~ 2 mK over characteristic veloc- 
ities of ~ 500 km s _1 . Similarly, the spectrum of W Comae 
shows flux variations of ~ 0.5 mK over characteristic velocities 
of ~ 500 km s _1 . For this reason, we did not attempt to model 
the baseline over the entire 4 GHz bandwidth, instead we focused 



our analysis on a velocity window of ±1000 km s _1 centered 
at the systemic redshift that is known from optical spectroscopy. 
In this interval, we subtracted a polynomial model fitted to the 
data from the spectrum. Then, we smoothed each spectrum with 
a Hann window and we binned the data to achieve a final resolu- 
tion of ~ 45 km s _1 per channel. Figure [2] shows the smoothed 
continuum-subtracted spectra for the three blazars. 

For RGB J07 10+591, we subtracted a linear baseline for both 
the 12 CO(l - 0) and the 12 CO(2 - 1) transitions. While fitting the 
baseline, we excluded a symmetric window of ±200 km s~ cen- 
tered at the systemic redshift to prevent the subtraction of weak 
emission lines. A higher order polynomial fit (order 10 for the 
12 CO(l - 0) line and order 5 for the 12 CO(2 - 1) line) was required 
instead for W Comae due to the lower quality of the baseline, again 
excluding a symmetric window of ±200 km s _1 around the op- 
tical redshift. We emphasize that, despite the high order adopted, 
the model for the baseline has characteristic shapes on scales of 
several hundred km s -1 and, therefore, it does not introduce spu- 
rious signals on small velocity scales (< 200 — 300 km s _1 ), 
preserving the noise properties of the original spectrum. Finally, 
for 1ES 1959+650, the 12 CO(2 - 1) emission line is visible at 
~ 100 km s _1 even before modeling the baseline. In this case, 
we adopted a polynomial fit of order 4 with a — 300 km s~ 1 win- 
dow. Conversely, the poor quality of the baseline at longer wave- 
lengths prevents us from detecting the associated 12 CO(l — 0) line 
trivially. After subtracting a baseline of order 15, a line is clearly 
detected when, based on the position of the 12 CO(2 — 1) line, a 
— 300 km s _1 window is excluded from the continuum fit (bot- 
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Figure 2. Spectra of the 12 CO(l — 0) and 12 CO(2 — 1) transitions, centered at the systemic velocity that is known from optical spectroscopy. The dashed and 
dash-dotted horizontal lines mark one and three times the rms of each spectrum. For 1ES 1959+650 (bottom panels), we show two different spectra that have 
been obtained with a fine-tuned window in the velocity interval — 300 km s — 1 (left) and a symmetric window ±200 km s — 1 (right). Gaussian models for 
the detected lines in 1ES 1959+650 are also shown with blue solid lines. 



torn left panel of Figure[2}. Instead, if we impose a symmetric win- 
dow of ±200 km s _1 around the optical redshift as previously done 
for the other two blazars, the 12 CO(l — 0) line is only marginally 
detected and the peak flux is suppressed by ~ 65% (bottom right 
panel of Figure [2). For this reason, we consider the 12 CO(l - 0) 
line flux affected by significant systematic uncertainty associated 
with the baseline subtraction. 

The root-mean-squared temperatures (T rms ) computed in the 
baseline-subtracted spectra over the ±1000 km s _1 velocity inter- 
val at the final resolution are listed in Table [T] in units of both the 
antenna temperature T% and of the flux density S v . For the flux 
calibration, we adopt the frequency dependent conversion S/T% 
available from the EMIR commissioning report. 



4 THE MOLECULAR CONTENT OF BL LAC OBJECTS 

Figure [2] shows that molecular gas is only detected in 1ES 
1959+650, i.e. in one out of three sources. Inspecting more closely 
the 12 CO(l - 0) spectrum for RGB J0710+591, a positive excess is 
visible at —300 km s _1 . Integrated over two channels, this spectral 
feature is significant at < 2a after accounting for both the statisti- 
cal error and the uncertainty of the baseline, which we estimate to 
be comparable to the measured T rms . For this reason, and because 
we do not see the corresponding 12 CO(2 — 1) transition, we only 
report an upper limit for RGB J0710+591. 

More significant, instead, is the detection for 1ES 1959+650. 
The 12 CO(2 — 1) line is detected at 3.2a accounting for both ran- 



dom and systematic errors, and it is not particularly sensitive to the 
baseline subtraction. Further, the 12 CO(l — 0) line is also detected 
at 4.4<r, but as previously discussed this transition is more sensitive 
to the baseline subtraction. Moreover, both lines appear at the same 
redshift (z co = 0.0475 ± 0.0002 and z co = 0.0474 ± 0.0002 for, 
respectively, 12 CO(l - 0) and 12 CO(2 - 1)) that is consistent with 
the one measured from optical lines. Finally, both line profiles are 
marginally resolved. 

For 1ES 1959+650, we model both lines with a Gaussian pro- 
file in the velocity interval 0—250 km s^ 1 . The peak of the molecu- 
lar emission is offset from the systemic velocity inferred from opti- 
cal spectroscopy by ~ 130 km s _1 , but since the uncertainty of the 
optical redshift is unknown, this shift may not be significant. The 
12 CO(1-0) and 12 CO(2- 1) lines have a full- width half-maximum 
velocity An = 121 km s _1 and Av = 110 km s _1 , respec- 
tively. These lines ar e narrow compared t o the typical line widths 
for luminous quasars (IScoville et al.l2003l) and more cl osely resem- 
ble th e typical velocity width of early-type galaxies dYoung et al.l 
1201 lh . By integrating over the Gaussian profile, we derive a line 
intensitjQ J co = 754 ± 170 mjy km s _1 for 12 CO(l - 0) and 
Jco = 1647 ± 511 mjy km s" 1 for 12 CO(2 - 1). Here, the er- 
ror is computed as SI co oc \/2NT rm sV ICB , where N is the num- 
ber of channels at resolution v las within the line width. We also 
attempt to estimate the uncertainty on the baseline by performing 



1 By directly integrating the data in the velocity interval — 250 km s 1 , 
we find similar values, consistent within the errors. 
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Figure 3 . (Lg/t) 12 CO(l — 0) line luminosity as a fun ction of redshift in: early type galaxies (orange crosses ; lYoung et al]|201 ll). massive galaxies (green 
triangles; Saintonge et al. 201 lj), quasar s (gray squares; IScoville et a7]|2003l), low luminosity quasars (gray X 's; Bertram et al. 2007), infrared-excess quasars 
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W Comae and RGB J0710+591. The right-end-side axis translates the 12 CO(l — 0) luminosity into the molecular hydrogen mass, assuming a = 4 Mq (K 
km s — 1 pc 2 )~ 1 . The black dotted and dashed lines represent the luminosity at constant flux density as a function of redshift for the IRAM 30m telescope and 
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multiple fits to the continuum with varying choices of polynomial 
orders and window positions. This simple exercise reveals that an 
additional factor \/2 should be included in the error budget. How- 
ever, we emphasize that this error estimate is still sensitive to the 
functional form adopted to model the ripples that are present in the 
baseline, and thus an additional systematic uncertainty may affect 
our final values. Only future interferometric observations will be 
able to quantify this error. 

Finally, the integrated intensity for the 12 CO(l — 0) line can 
be converted into the luminosity L' co (in K km s" 1 pc 2 ) by using 

Ico ( D L \ 2 1 



l'co = 2.45 x irr 



Jy km s" 1 \Mpcy l + z 



(1) 



Assuming a typical conversion factor from CO luminosity to 
molecular hydrogen mass a = 4M (K km s" 1 pc 2 )" 1 (see 
lEvans et al.l200lh . we infer a total hydrogen mass Mh 2 = (3.04 ± 
0.69) xlO 8 M .A summary of the molecular gas properties of 1ES 
1959+650, together with the la upper limits for RGB J0710+591 
and W Comae, is provided in TableQ] 



5 DISCUSSION 

With two tight upper limits and one detection at Mh 2 ~ 3 x 
10 8 Mq, these pilot observations suggest that VHE-detected BL 
Lac objects do not generally bear large amounts of molecular gas. 
However, the detection of weak lines in 1ES 1959+650 proves that, 
despite the bright synchrotron continuum, deep enough observa- 
tions of the molecular lines can be a viable technique to constrain 
and refine the redshifts of sources of particular interest for which 



approximate distances are know from indirect methods (see Section 



5.1 A comparison with other AGN and galaxy populations 

Although our sample is too small for any conclusive analysis, in 
the left panel of Figure|3]we attempt a first comparison between the 
CO luminosity of BL Lac objects and other sources for which mil- 
limeter observations are available from the literature. All the data 
have been homogenized to a single conversion factor a — 4Mq 
(K km s _1 pc 2 ) -1 and all the limits are reported at 3<r, preserving 
the choice of the velocity width as published in the original studies. 
For the upper limits of RGB J07 10+591 and W Comae, we adopt a 
line width of 150 km s _1 that is motivated by the detection in 1ES 
1959+650. 

The difference between VHE-detected BL Lac objects and 
quasars dScoville et alj 



galaxies (Bertr am et al 



2003; lEvans et all 1200 lh or Seyfert I host 



20071) is evident, given that the molecular 



gas detected in quasars is up to a factor of ten higher than what was 
found in blazars by our observations. This discrepancy is appar- 
ent even if we consider only quasars from the lScoville etai1 ( l2003h 
sample that overlap in redshift with our blazars and that were se- 
lected purely based on optical luminosity. In fact, there is an even 
more pronounced separation between blazars and infrared ultra- 
lumi nous quasars tha t have a median molecular gas mass of 6 x 10 9 
M dXia et al^OlU) . 

Further, BL Lac objects are only consistent with the lower- 
end of the molecular gas distribution in galaxies that are selected 
purely based on their stellar mass (M star > 10 10 Mp)) regardless 
to the presence of AGN and morphological type dSaintonge et al.l 



© xxxx RAS, MNRAS 000,[T]i] 



6 Fumagalli et al. 



Conversely, the CO luminosity in BL Lac objects seems con- 
sistent with the level of molecular gas typically detected in radio- 
loud galaxies. These sources are indee d observed to have a lower 
molecular gas content th an quasars (Ocana Flaquer et all |2010| : 
ISmolcic & Riechersll20Tll) . Finally, VHE-detected BL Lac objects 
also appear comp atible with the population of early type galaxies 
dYoungetal]|201lh in which only a fifth of the observed sources 
are detected at the current sensitivity limits. It should be noted that 
these classes of objects are observed at lower redshifts than the 
blazars and the luminous quasars, thus the Malmquist bias com- 
plicates a systematic comparison between different classes. For in- 
stance, the subset of radio galaxies with z > 0.05 exhibits a molec- 
ular mass that is systematically higher and comparable to quasars 
than in the f ull sample (see ad ditional d iscussion on this redshift de - 
pendence in I Ocana Flaquer e t al. 2010; ISmolcic & Riechersl2011f) . 
Nevertheless, the comparison between quasars and BL Lac ob- 
jects is done at comparable redshifts. Further, because blazars lie 
at z > 0.05, we can still conclude that blazars are consistent with 
the median molecular mass content of early type galaxies and ra- 
dio galaxies given our current sensitivity limit. Deeper observa- 
tions will investigate whether this is true or whether blazars ex- 
hibit in fact an intermediate level of molecular gas between those 
of quasars and early type galaxies. 

The qualitative trends discussed above are confirmed by the 
more quantitative analysis presented in the right panel of Figure 
[3] where for each class we compare the 25 th , 50 th , and 75 th 
percentiles of the observed L'co distributions. For this , we use 
the Kaplan-Meier estimator ( Feigels on & Nelsor] 1 19851) for ran- 
domly censored dat a as implemented in the ASURV package 
dLavallev etaljl992h . Since our sample includes only three blazars, 
percentiles cannot be defined and we superimpose the observed 
data points in Figure[3] We can still derive a mean value of L' co = 
7.6 x 10 7 K km s _1 pc 2 that, however, is not particularly infor- 
mative due to small statistics. We further emphasize that the per- 
centiles shown in the right panel of Figure [3] refer to samples that 
are selected to be representative of the populations in exam, al- 
though they are not complete. 



5.2 Possible observational biases 

Noting the lack of strong emission from these BL Lac objects, it 
is worth examining whether the discrepancy between this class of 
blazars and quasars can be attributed to an observational bias and, 
in particular, whether the bright synchrotron continuum may out- 
shine strong emission lines. 

To rule out this possibility, we constructed mock spectra with 
different continuum intensities and line properties similar to the one 
measured in the 12 CO(l - 0) transition for 1ES 1959+650. Follow- 
ing the same procedure adopted for real data, we analyzed these 
mock spectra multiple times, adding different realizations of Gaus- 
sian noise to a level that is comparable to the observed one. At the 
end of this test, we were able to recover the input line parameters 
with a dispersion that is consistent with the estimated uncertainties. 
We also built an additional set of mock spectra, with the noise prop- 
erties inferred for RGB J0710+591. Superimposing to the contin- 
uum a line of total intensity equal to three times T rms , we were able 
to visually identify the emission line. We conclude that, if present, 
emission lines significantly brighter than the quoted limits should 
have been detected despite the bright synchrotron continuum. This 
procedure does not take into account the ripples that are present at 
3 mm in the spectra of W Comae and 1ES 1959+650, and future 



interferometric observations will assess the error introduced by the 
poor quality of the baseline. 

Next, we consider whether the different luminosities in BL 
Lac objects and quasars can be attributed to projection effects. 
Since the emitted flux from optically thick regions is proportional 
to the gas surface area, different viewing angles of the same source 
can result in different L' co . Assuming that the gas settles in a disk 
and that the jets extend above the plane, observations of blazars 
probe the face-on view of the emitting disk. Thus, in order to 
explain the lower luminosity of BL Lac objects, the gas should 
be concentrated in a small region where the disk surface area is 
comparable or smaller than the disk scale height. However, this 
seems implausible for at least two reasons. First, a concentration 
of 10 8 M© of molecular gas in the inner few hundred parsecs 
close to the black hole would produce a hydrogen density that 
exceeds what is inferred fro m the narrow and broad line regions 
dOsterbrock & Fe rland 2006), by at least two orders of magnitude. 
Second, the few observations of quasars where the gas disk is re- 
solved show that the emitting regi ons typically extend to ~ 1 kpc 
or, in some cases, up to a few kpc dRiechers et alj|201 11) . 

5.3 On the difference between quasars and BL Lac objects 

Cons i dering morpholo gical studies of blazars (e.g. IScarpa et al.l 
l2000l : lUrrv et al1l2000h . a lack of molecular gas is not completely 
surprising since almost all the BL Lac objects resides in massive 
early- type galaxies. These studies further suggest that the nuclear 
activity is not necessarily triggered by recent interactions that have 
a dramatic effect on the large-scale morphological and structural 
properties of the host galaxies. In fact, radio galaxies, blazars and 
radio-quiet elliptical galaxies share similar structural and morpho- 
logical characteristics. Our CO observations add additional evi- 
dence in support to this idea: in our limited sample the cold gas 
in BL Lac objects is consistent with the lower level of CO seen in 
early- type and radio galaxies (see Figure [3). 

Further, low-excitation radio AGN (LERAGN) and Fanaroff 
and Riley class I (FR-I) sources have a lower molecular gas con- 
tent than high-excitation radio AGN (HERAGN) and FR -II sources 
dSmolcic & Riecherdl2m~il ; lOcafia Flaquer et al.ll2010h . Although 
it is still debated whether the discr epancy between FR-I and F R-II 
sources is due to selection effects (Ocana Flaque ret alfeOloT) . the 
low molecular content found in BL Lac objects is consistent with 
these findings, if we allow for an admittedly loose identification of 
LERAGN with FR-I class, which contains most of the BL Lac pop- 
ulation. Thus, based on the properties of their interstellar medium, 
low-redshift quas ar and BL Lac host g alaxies appear as two distinct 
families (see also lScoville et al . 2003) that are not drawn from the 
same parent population. 

In light of recent findings that both th e amount of atomic inter- 
stellar medium (ISMi lFabello et alJlioTlh an d the bulk velocity of 
the molecular content in nearby AGN hosts Jciuillard et al.l l2012h 
are only weakly affected by the nuclear activity, the low molecular 
content in BL Lac objects might be a consequence of the quenching 
mechanisms associated with the growth of bulges an d an increase 
in the stellar surface density dKauffmann et al.l 20121) rather than a 
direct consequence of the presence of a jet. In any case, regardless 
of whether the AG N are causally linked to the gas properties of the 
host galaxies (e.g. iNesvadba et alj|2~01 ll) . the ISM content of BL 
Lac objects revealed by our observations is typical of late stages of 
galaxy evol ution, as it is also advocated to explain the blaz ar se- 
quence (e.g. iGhisellini et alj[l998l : iBottcher & Dermerll2002l) : BL 
Lac objects correspond to late stages of activity in which the cen- 
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tral black hole is fed by a radiatively inefficient accretion disk. This 
view is further supported by the fact that VHE emitting HSP BL 
Lacs are generally found to be best represented by simple SSC 
models, while VHE emitting ISP and LSP objects require external 
photon fields to describe their broadband data. Moreover, mid-IR 
observations of a large population of BL Lac objects show no clear 
signatures of torus emission, con sistent with the expected low ac- 
cretion rates dPlofkin et al.ll2012h . Thus, our observations fit within 
unified models for AGN in which the accretion rate shapes the ob- 
served nuclear properties. BL Lac objects are not simply AGN with 
radio emission superimposed, as described in unification model s 
that are based solely on geometrical obscuration dAntonuccill 19931) . 
but are objects in which the onset of the radiatively inefficient ac- 
cretion flow has led to the disappe arance of the broad emission lines 
fe.g. lHcl|200&|Trump et ai1l201lh . 

In closing, we note that the discrepancy in CO luminosity be- 
tween BL Lac and quasar host galaxies suggested by our obser- 
vations has consequences that extend beyond the study of AGN. 
It has recently been noted that plasma instabilities that develop in 
the jets of TeV-emitting blazars could inject a considerable amount 
of heat in the intergalactic medium {Broderick e t al.ll20lH) . alter- 
ing the therm al history of the Universe at low redshifts and at 
low densities dChang etalj|201lb . However, if BL Lac objects and 
quasars are not drawn from the same parent population, the num- 
ber density of blazars and quasars may not evolve in parallel across 
all redshifts. Thus, the underlying assumption that the TeV blazar 
activity does not lag that of quasars may need additional verifi- 
cation. Additionally, observations of high-excitation radio galax- 
ies show a stronger cosmic evolution at all luminosities compared 
to th e low-excitation radio galaxies which e xhibit milder evolu- 
tion dWillott et alj l200ll; ISadler et al.l |2007|; iDonoso et all 120091 ; 
ISmolcic et al.l 120091 : iBest & Heckmanl l2012h . A similar trend is 
also e xpected for the relative evolution of FSRQs and BL Lac ob- 
jects dCavaliere & D'EliabOO^lBottcher & Dermerl2002l) . Besides 
this specific application, understanding the connection between the 
heating mechanisms of blazars and the ISM properties of the low- 
excitation radio galaxies becomes relevant in the context of cur- 
rent feedback prescriptions such as the "radio-mode" f eedback that 
is commonly adopted in semi-analytic models (e.g. ICroton et al.l 
2006: lBoweretalJl2006l) . 



6 SUMMARY AND FUTURE PROSPECTS 

BL Lac blazars detected at VHE are key sources to probe a variety 
of astrophysical processes, including the physics of jets, the prop- 
erties of the EBL and the strength of the IGME However, due to 
their featureless optical spectra, it is often impossible to establish a 
redshift with optical emission or absorption lines, limiting the use 
of these objects in cosmological studies. Motivated by the need to 
reliably determine distances for VHE BL Lac objects with other 
techniques, we have undertaken a pilot program at the IRAM 30m 
telescope to assess the molecular gas content in three sources with 
known redshifts to establish whether molecular lines can be used to 
trace distances for this class of blazars. 

For one BL Lac object (1ES 1959+650), we detected both the 
12 CO(l - 0) and the 12 CO(2 - 1) transitions at > 3cr confidence 
level, corresponding to a molecular mass of Mh 2 ~ 3 x 10 8 Mq. 
However, due to the poor quality of the baseline, we regard this 
detection as affected by a large systematic uncertainty. Conversely, 
for the other two sources (W Comae and RGB J0710+591), we do 
not detect CO emission to 3<j limits L' co < 2 x 10 8 K km s _1 pc 2 



and L' co < 4 x 10 8 K km s _1 pc 2 , respectively. In both cases, we 
assumed a line width of 150 km s _1 , motivated by the lines de- 
tected in 1ES 1959+650. Assuming a conversion factor a = 4Mq 
(K km s _1 pc 2 ) -1 , these limits translate to Mh 2 < 8 x 10 8 M Q and 
Mh 2 < 1.6 x 10 9 M Q . 

Although this sample is too small to derive robust statistical 
conclusions, our observations suggest that BL Lac host galaxies 
may not, in general, bear large amounts of molecular gas. The de- 
tected levels are in fact comparable to what is observed in early- 
type or radio galaxies. Thus, a search for molecular lines is not a 
feasible approach to establish blazar redshifts with instruments that 
offer a limited bandwidth. Nevertheless, current facilities such as 
the IRAM Plateau de Bure Interferometer or the Atacama Large 
Millimeter Array can reach up to one order of magnitude better 
sensitivity than what is achievable in similar integration times with 
the IRAM 30m telescope (see Figure[3]l. Therefore, the detection of 
molecular lines can still be a promising way to obtain precise red- 
shifts for blazars whose approximate distances have already been 
estimated through indirect methods such as the ab sorption from the 
intervening Lya forest (e.g.lDanforth et al.l l2010). the EBL attenu- 



ation (e.g. Steck er & Scullvl2010h, or th e photometric properties of 



host galaxies (e.g. lSbarufatti et al .120051) . 

Besides the possibility of establishing redshifts, these obser- 
vations offer tantalizing hints that the molecular gas content of BL 
Lac objects is lower than that of quasars. Larger samples are now 
needed to assess the significance of this discrepancy, but a sys- 
tematic deficiency of molecular gas in BL Lac objects compared 
to quasars will provide confirmation that the relative host galax- 
ies are not drawn from the same parent population. Our observa- 
tions, together with previous studies of blazars and radio sources 
(e.g. the low and high excitation radio AGN or the FR-I and FR-II 
sources), provide increasing evidence that BL Lac objects represent 
later stages of galaxy evolution, in which the central black hole is 
fed by a radiatively inefficient accretion disk. In the era of large 
millimeter arrays, a multiwavelength view of the ISM of AGN is 
required to fully explore the connection between different types of 
active nuclei and the host galaxies. 
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